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Abstract 

We give a notion of BV function on an oriented manifold where a volume form and 
a family of lower scmicontinuous quadratic forms G p : T p M — > [0,oo] are given. When 
we consider sub-Riemannian manifolds, our definition coincides with the one given in 
the more general context of metric measure spaces which are doubling and support a 
Poincare inequality. We focus on finite perimeter sets, i.e., sets whose characteristic 
function is BV, in sub-Riemannian manifolds. Under an assumption on the nilpotent 
approximation, we prove a blowup theorem, generalizing the one obtained for step-2 
Car not groups in ptjFI 
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1 Introduction 

Sub-Riemannian manifolds are a class of length spaces of non-Euclidean type having a 
differentiable structure. Our interest in studying functions of bounded variation in this 
framework arises from the aim of understanding the structure of finite perimeter sets in 
the general sub-Riemannian setting. This clearly requires suitable notions of "intrinsically 
regular" hypersurfaces, rectifiability, reduced boundary and blowups. 

Sobolev and BV functions have been investigated in W 1 endowed with the Lebesgue 
measure and with the Carnot-Caratheodory distance d cc associated with a family of vec- 
tor fields. Under the assumption that the family is Lie bracket generating, the Lebesgue 
measure is doubling with respect to the Carnot-Caratheodory distance [37], a Poincare 
inequality holds [30] and (M n ,d cc ) is complete [16]. These are the main assumptions which 
enable the authors in [25] to establish Sobolev and isoperimetric inequalities as well as an 
approximation theorem of Meyers-Serrin type (see also [22] for a related result with weaker 
regularity assumptions on the vector fields). 

Our main goal is to develop a systematic theory of BV functions and sets of finite 
perimeter in manifolds with suitable structures. To this aim, one needs two ingredients: 
first a volume measure (with respect to which an integration by parts formula will hold); 
second a notion of length of tangent vectors (along which one calculates distributional 
derivatives). For the volume measure, when a manifold M is oriented, it suffices to take 
a non-degenerate n-form uj which induces the orientation of M (where n = dimM) and 
consider the measure m defined on Borel sets B C M by m(B) = J B uj. Given an open set 
n C M and a vector field X, a function u 6 L 1 (17,m) has distributional derivative along 
X if there exists a Radon measure Dxu on £1 such that 



where div^X is defined in ([5]). To define the length of tangent vectors, we use a family of 
lower semicontinuous quadratic forms G p : T p M — > [0, oo] defined on the tangent bundle 
of M. Note that the dimension of the vector space T>(p) = {v € T p M \ G p (v) < oo} may 
vary with respect to the point. Taking this into account, it is natural to say that a function 
u £ L 1 (0,m) has bounded variation in f2 if, for every smooth vector field X such that 
G p (X(p)) < 1, p £ f2, the distributional derivative Dxu is a Radon measure of finite total 
variation in 0, and 



the supremum being taken among all smooth vector fields such that G(X) < 1 on £1. 
Thus, we write u £ BV(£l, g,uj), where g is the scalar product associated with G. More 






D g u\\(Q,) := sup \D x u\(Cl) < oo 



(1) 
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precisely, for each p G M we have that g p (-,-) is the unique scalar product on V{p) such 
that g p (v,v) = G p (v) for every v G D(p), Section [2.31 for more information. 

In this quite general setting, distributional derivatives can be weakly approximated by 
derivatives of smooth functions, that is, a Meyers-Serrin theorem holds (see Theorem 12. 4p . 
Moreover, the fact that D^u is a Radon measure with finite bounded variation is charac- 
terized in terms of difference quotients along the flow generated by X (see Theorem 12.51) . 
Using G, one can define the Carnot-Caratheodory distance d cc between points of M as 
the infimum of lengths of absolutely continuous curves connecting the two points, where 
length of tangent vectors is computed with respect to G. When d cc is finite, (M, d cc , m) 
is a metric measure space. It is then natural to compare the space BV(Q,g, oj) with the 
notion of BV function in a metric measure space developed in [36] (see also [6]). Without 
further assumptions, W6 CcLii only show that BV(^l^ d cc ^ m) is embedded in BV(p.,g,uS) (see 
Theorem 12. 7p . 

Oriented sub-Riemannian manifolds, where G is induced locally by bracket generating 
families of vector fields, cast in the framework above. In this case, on coordinate charts, G 
is given by 



(with the convention inf0 = oo) where X±, . . . , X m play the role of orthonormal frame 
and dim(Xi(p), . . . ,X m (p)) may vary with respect to p. In particular, the aforementioned 
notion of BV space encompasses the classical one in oriented Riemannian manifolds, the one 
associated with a Lie bracket generating family of vector fields in M. n and it also includes the 
rank- varying case, e.g. the Grushin case and almost-Riemannian manifolds (see Section [3TTj) . 

In this setting, the approximation result (Theorem 12. 4p allows to show that the met- 
ric and differential notion of bounded variation coincide and the corresponding spaces 
BV(£l, g, uj) and BV(£l, d cc , m) are isometric (Theorem 13. ip . A first consequence of this fact 
is that the set function ||-D g u|| defined on open sets as in ([T|) is a Borel measure. Moreover, 
the Riesz theorem of Euclidean setting (see for instance [201 Theorem 1 page 167]) can be 
generalized. More precisely, if u G BV(£l, g, u>), then there exists a Borel vector field v u 
satisfying G{v u ) = 1 ||Z) fl it||-a.e. in VL. Moreover, for every smooth vector field X such that 
G(X) < 1 in 0, the distributional derivative of u along X can be represented as 



Without assumptions on the dimension of T>(p), even if a local basis X±, . . . ,X m inducing 
G is given, some care is needed, due to the fact that a smooth vector field X satisfying 
G(X) < 1 is in general a linear combination of the Xi with coefficient in L°° only. 

When we consider sets of finite perimeter, i.e., sets whose characteristic function has 
bounded variation, this result provides a notion of geometric normal (which corresponds to 
its Euclidean analogue) and which is a horizontal Borel vector field of G- length 1. 

In Euclidean metric spaces, the structure of finite perimeter sets has been completely 
understood since De Giorgi's seminal works |18[ll9j. In this context, if E has finite perime- 
ter, then the perimeter measure of E is concentrated on a set which is rectifiable and it 




Dxu = g(X,u u )\\D g u 



(2) 



3 



has codimension one. The main step behind this result is a blowup theorem, showing that 
when p belongs to the reduced boundary of E, the sequence of blowups E r = 5\/ r {E + p) 
converges to a halfspace in L\ . 

In non-Euclidean metric spaces, after [31], rectifiability theory has been developed in 
Banach spaces [T2J Q3], and in Carnot groups [231 E51 EB] , which are Lie groups whose algebra 
admits a stratification with respect to a one parameter group of dilations 6 r . More precisely, 
in [24] the authors generalize De Giorgi's theorem in Carnot groups of step 2. Their proof is 
inspired by the one in the Euclidean case. Moreover, Carnot groups are homogeneous and 
this makes the perimeter measure both homogeneous with respect to dilations and invariant 
by translations. Joining these properties with isoperimetric inequalities along with the 
compact embedding in BV shows that bounded sequences of rescaled sets have converging 
subsequences and the blowups are both monotone along a horizontal direction and invariant 
along all orthogonal directions. The techniques of [23] have been further extended in [T7] 
to a special class of Lie groups that do not possess dilations, see Example 13. 3[ where the 
"linearization" of the left invariant vector fields is obtained by the Rothschild-Stein lifting 
theorem, |39j . 

As a first step toward a generalization of De Giorgi's theorem in sub-Riemannian mani- 
folds, in this paper we show a blowup theorem which generalizes the one |24|, Theorem 3.1] 
in Carnot groups of step 2. Again, the proof of Theorem 14. 21 is inspired by the corresponding 
one in the Euclidean case. However, the rationale behind our proof is somehow different 
from the one in [24] , Without a Carnot group structure, the main idea is to exploit two well 
known facts in sub-Riemannian geometry |15] : a metric tangent cone to the manifold at a 
point p always exists; the quasi-isometry between dilated balls centered at p and balls in 
the metric tangent cone can be given explicitly by a system of suitable coordinates (called 
privileged, see Definition 13. 3D centered at p and, in particular, it is a diffeomorphism ip p . In 
this coordinate system, there exists a subgroup of dilations 8 r intrinsically associated with 
the sub-Riemannian structure (and centered at p). Hence, given a finite perimeter set E 
and p in its reduced boundary (see Definition I3.2[) . reading E through ip p , it makes sense 
to consider the blowups E r = 8i/ r <p p (E). Our result states that if the metric tangent cone 
to (M, d cc ) at p is a Carnot group of step 2 then 

LjJ^-limlE!,. = If, 

where F is the vertical halfspace in the Carnot group associated with the geometric normal 
ve(p) (for the precise statement, see Theorem I4.2L ) In particular, we are able to show that 
the sequence of blowups {l_E r } r >o is compact in L\ oc and that if lp is a L^-limit then F 
is monotone along the geometric normal ve{p) and invariant along orthogonal directions 
to ve{p)i in the distributional sense. To prove compactness, we exploit the fact that the 
distance in the metric tangent cone is the limit of Carnot-Caratheodory distances associ- 
ated with a sequence of "perturbed" vector fields (see Theorem I3.5p . Properties of limits 
are consequences of the definition of geometric normal and of the asymptotically doubling 
property of the perimeter measure. Finally, it is only at this stage of the proof that we 
invoke the fact that the metric tangent cone at p is a Carnot group of step 2, to show that 
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the limit of the rescaled sets actually exists. The latter assumption is fulfilled of course 
when the manifold itself is a Carnot group of step 2 but also in more general case, e.g. in 
the step 2 equiregular case (see also Example 13 .7|) . This hypothesis is essential as it was 
shown in [23] that in Carnot groups of step higher than 2 the blowup at a point in the 
reduced boundary need not be a vertical hyperplane. Without bounds on the step, the only 
result available so far is [13], where it has been proved that at almost every point (with 
respect to the perimeter measure) there exists a subsequence of blowups converging to a 
vertical halfspace. 

Let us mention an application of our result in the rank-varying case. As we see in 
Example 13.51 there exist sub-Riemannian manifolds with the property: for every point p 
the metric tangent cone at p is either the Euclidean space or a Carnot group of step 2. 
For these manifolds, which are also called almost-Riemannian, the horizontal distribution is 
rank-varying and it has full rank at points where the tangent cone is the Euclidean space. 
Combining our theorem with the one in the Euclidean case, we obtain that, for sets of finite 
perimeter in these manifolds, the blowups at points in the reduced boundary converge to a 
halfspace. 

Another important corollary of our blowup theorem is that, setting h{B r (p)) = ITE^rM) 
where B r {p) is the open ball for d cc , the density 

r \\D g ±E\\{B r {p)) 

11111 Ul-D I \\ ( 3 ) 

exists and equals to the perimeter of F in the unit ball in the Carnot group divided by the 

l E ||(B r (p)) 



Lebesgue measure of the unit ball. This improves the weaker estimates on - — i h ^ B ^ — 
which have been proved in Theorem 5.4] in the metric setting. Moreover, denoting by 
S h the spherical measure build by the Caratheodory's construction (see |2H 2.10.1]) with 
h as gauge function, the existence of the limit in ([3]) implies upper and lower bounds on 
the Radon-Nikodym of ||£) fl ls|| with respect to S h (restricted to the reduced boundary of 
E). As the Radon-Nikodym derivative of the perimeter measure with respect to 

S h has been shown to exist in the metric context (see [3 Theorem 5.3]), an open question 
is whether this derivative coincides with ([3j) for ||Z) fl l^||-almost every p. In the constant 
rank (equiregular) case, a related result in [2] computes the density of the spherical top- 
dimensional Hausdorff measure (where Q is the Hausdorff dimension of any ball) with 
respect to m in terms of the Lebesgue measure of the unit ball in the metric tangent cone. 

The paper is organized as follows. We define distributional derivatives along vector fields 
in manifolds with a volume form in Section 12.21 Using a family of metrics in the tangent 
bundle we then define the space of BV functions in Section 12.31 and we prove an approxi- 
mation results for distributional derivatives. In Section [2.41 we show that BV{Sl, d cc , m) is 
continuously embedded in BV(£l, g,u). Section [3] is a primer in sub-Riemannian geometry. 
Even though our main results are local, we find it useful to recall the general definition of 
sub-Riemannian structure that relies on images of Euclidean vector bundles and includes 
the rank- varying case. In Section [3.11 we list some significative examples, including Carnot 
groups. In Section 13.21 we analyze the notion of BV space in sub-Riemannian manifolds. 
First, we specify its relation with the BV space defined in a metric measure space, showing 



5 



that the two BV spaces are actually isometric. Second, we prove a Riesz theorem which 
generalize the Euclidean analogue. Section 13.31 recalls the notion (and basic properties) of 
privileged coordinates and nilpotent approximation (for more details we refer the reader to 
|15j.) In Section [3.41 we explain the relation between nilpotent approximations and metric 
tangent cones to sub-Riemannian manifold at a point and, under an additional assump- 
tion at the point, we recall how to show that the nilpotent approximation is isometric to a 
Carnot group endowed with the control distance induced by a left invariant metric on the 
horizontal bundle. In Section H] we prove the blowup theorem. We split the proof into two 
main parts. In Sections 14.11 and 14.21 we demonstrate compactness of the dilated sets and 
monotone and invariance properties of limits. Then, in Section [4.31 we use the assumption 
on the nilpotent approximation to provide the existence and characterize the limit of dilated 
sets vertical halfspace. 

2 Preliminaries 

2.1 Basic notation and notions 

Given a set A C B, we will use the notation 1^ : B — > {0, 1} for the characteristic function 
of A, equal to 1 on A and equal to on B \ A. In a metric space (X, d), the notation B r (x) 
will be used to denote the open ball with radius r and centre x. 

Differential notions. Throughout this paper, M denotes a smooth, oriented, connected 
and n-dimensional manifold, with tangent bundle TM. The fiber T X M can be read as the 
space of derivations on germs of C 1 functions cp at x, namely [f</9](x) = dip x (v), for v £ T X M. 
In the same spirit, we read the action of the differential dF x : T X M — > T X N of a C 1 map 
F : M ->• N as follows: 

df x (v)(tp) =v{(poF){x) ipeC^N). 

Any C 1 diffeomorphism F : M — > N between smooth manifolds induces an operator 
: TM TN, by the formula = dF x (X(x)); equivalently, in terms of 

derivations, we write 

[(F*X)<p] o F = [Xfa o F)] V<peC l (N). 

Measure-theoretic notions. If J 7 is a u-algebra of subsets of X and \i : T — > W 71 is a 
<T-additive measure, we shall denote by \fj,\ : T — > [0, oo) its total variation, still a u-additive 
measure. By the Radon-Nikodym theorem, fi is representable in the form g\n\ for some 
^-measurable function g : X — > M. m satisfying \g(x)\ = 1 for |/i|-a.e. x £ X. Given a Borel 
map F, we shall use the notation F# for the induced push-forward operator between Borel 
measures, namely 

F # fi(B) := /j,(F~ 1 (B)) for all B Borel. 
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2.2 Volume form, divergence and distributional derivatives 

We assume throughout this paper that M is endowed with a smooth re-form u). We assume 
also that J M foj > whenever / G C\{M) is nonnegative and not identically 0, so that the 
volume measure 

m{E) = u, EC M Borel (4) 
Je 

is well defined and, in local coordinates, has a smooth and positive density with respect to 
Lebesgue measure. Accordingly, we shall also call uj volume form. 

The volume form uj allows to define the divergence div^X of a smooth vector field X as 
the smooth function satisfying 

diw^X uj = L x u, (5) 

where Lx denotes the Lie derivative along X. Using properties of exterior derivative and 
differential forms, we remark that div^X is characterized by 

- f <pd\v ul Xuj = [ (Xip)u Vtp£Cl(M). (6) 
JM JM 

By applying this identity to a product ftp with / G C 1 (M) and tp 6 C£°(M), the Leibnitz 
rule gives 

-/ ftpdiVuXco- f f(X<p)uj= [ tp(Xf)u V^GC C 00 (M). (7) 

JM JM JM 

We can now use this identity to define Xf also as a distribution on M, namely Dxf = g 
in the sense of distributions in an open set C M means 

- f f^div^Xu- [ f{Xy)u= [ ipgu \/ip G C£°(f2). (8) 
Jn Jn Jn 

Our main interest focuses on the theory of BV functions along vector fields. For this 
reason, we say that a measure with finite total variation in f2, that we shall denote by Dxf, 
represents in f2 the derivative of / along X in the sense of distributions if 

- f fipdw^Xu- [ f(X(p)u= [ cpdDxf V99GC c °°(fi). (9) 
Jn Jn Jn 

By g]) and ©, when / is C 1 we have D x f = (Xf)m. 

We can now state a simple criterion for the existence of Dxf, a direct consequence of 
Riesz representation theorem of the dual of C c (0). In order to state our integrations by 
parts formulas ([8]), ([9]) in a more compact form we also use the identity 

div w ((pX) = (pdiv^X + Xip. 

Proposition 2.1. Let £2 C M be an open set and f £ L} (Cl,m). Then Dxf is a signed 
measure with finite total variation in £1 if and only if 

sup ij f d\v w ((pX) uj \<peC™(n), |^|<l|<oo. 

If this happens, the supremum above equals \Dxf 
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A direct consequence of this proposition is the lower semicontinuity in Lj oc (ft, m) of 
/ I—?- \Dx f\(Q)- We also emphasize that, thanks to ©, we have the identity 

D^ x f = i;D x f V*/>GC°°(ft), (10) 

and the properties ([8]) and © need only to be checked for test functions if whose support 
is contained in a chart. 

2.3 Distributions, metrics and BV functions on manifolds 

On M we shall consider a family of lower semicontinuous quadratic (i.e. 2-homogeneous, 
null in and satisfying the parallelogram identity) forms G x : T X M — > [0, oo] and the 
induced family V of subspaces 

V(x) := {v G T X M | G x (v) < oo} . 

We are not making at this stage any assumption on the dimension of T>(x) (which need 
not be locally constant) or on the regularity of x t— > G x . We shall only assume that the 
map (x,v) \— > G x {v) is Borel. This notion can be easily introduced, for instance using local 
coordinates. Obviously G x induces a scalar product g x on V(x), namely the unique bilinear 
form on T>{x) satisfying 

9x{v,v) = G x (v) VveV(x). 
For SJcM open, we shall denote by r(ft,P) the smooth sections of T>, namely: 

T(n,V) :={X\X is smooth vector field in ft, X(x) G V(x) for all iGS!}. 

We shall also denote 

T 9 (n,V) := {X G T(Q,V) | g x (X(x),X(x)) < 1 Vx G 0} . 

Note that both T> and g are somehow encoded by G. Thus the following definition of 
BV space only depends on G and u. 

Definition 2.1 (Space BV(Q, g,uj) and sets of finite perimeter). Let ft C M be an open set 
and u G L 1 (ft,m). We say that u has bounded variation in ft, and write u G BV(£l,g,uj), 
if D x u exists for all X G Tf (ft, V) and 

sup{|£>xH(ft) | X G r 9 (ft,£>)} < oo. (11) 

We denote by ||-Dpu|| the associated Radon measure on ft. If E C M is a Borel set, we say 
that E has finite perimeter in ft if 1^; G BV(£l, g, u>). 

Remark 1. Let us point out that replacing ft in (lllj) with any of it open subsets, we get 
a set function on open sets. Thus, the so-called De Giorgi-Letta criterion, see for instance 
|1H Theorem 1.53], allows us to extend this set function to a Radon measure on ft. 
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Equivalently, thanks to Proposition 12. 1\ we can write condition (jlip as follows: 

supjy udiv u (<pX)u | X G r»(fi,X>), € C~(M), |p| < l| < oo. (12) 

These .BV classes can be read in local coordinates thanks to the following result (more 
generally, one could consider smooth maps between manifolds). 

Proposition 2.2. Let Q C M and U C W 1 be open sets, let u G L 1 ^) and let <f> G C°°(Cl, U) 
be an orientation-preserving diffeomorphism with inverse tp. For each y G U, set 

G y {w) := G m (d y iP(w)) yw G T X U ~ M. n , u(y) := u^(y)) (13) 

and define T> in U and a metric g inT> accordingly. Then, setting uj := i^*uj, the following 
holds 

<f>#( D xu) = D^ x u vier(fi,D) (14) 
in the sense of distributions. In particular u G BV(Q, g, oj) if and only if u G BV(U,g,u>). 

Proof. Let X be a smooth section of T> with compact support contained in SI and consider 
the change of variable 

/ u div^X oj = / u (div^X) o tp oj . 
Jn Ju 

By Lemma 12.31 below, it follows that 



u div^X oj = u (Hvq(4>*X) oj, 
Q Ju 

where we denote by div^y the divergence of Y with respect to oj. If we apply this identity 
to the vector field ipX and use the relation 4>*((pX) = ((p o (fi~ 1 )cj)*X we obtain (JT3|) . 
Finally, since for all y G U we have G y (4>*X(y)) < 1 if and only if 

G m {d y ^X(y))) = G m {X{^{y))) < 1 

the claim follows. □ 

Lemma 2.3. Under the assumptions of Proposition PO?l the following formula holds 

(div w X)o^ = div (^JT). (15) 

Proof. If ip G C£°(fJ), then a change of variable in the oriented integral yields 

ip div^X uj = Xtp to = (Xip) o ij) oj = / [(^X)^] o ip oj , 
n Jn Ju Ju 

where the last equality is a consequence of the definition X = <f>*X. The previous equalities 
give 

— ip div^X uj = I X((p o ip) oj = — / (p o ip div^X oj = — ip (div^X) o <f> oj . 
Ju Ju Ju Jn 
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The arbitrary choice of <p proves the validity of (|15p . □ 

Distributional derivatives of L 1 functions as defined in (J9]) can be weakly approximated 
by derivatives of smooth functions as we prove in the next theorem. 

Theorem 2.4 (Meyers-Serrin). Let Q C M be open, u G L 1 (Q,m). Then there exist u n G 
C°°(Q) convergent to u in Z/j (12, 7TJ») satisfying Xu n 

weakly, in the duality with C C (Q), for all smooth vector fields X in £1 such that Dxu exists. 
The same is true if we consider vector-valued measures built with finitely many vector fields. 

Proof. By a partition of unity it is not restrictive to assume that O is well contained in 
a local chart. Then, by Proposition 12.21 possibly replacing g and uj by their counterparts 
g and uj, we can assume that C W n . In this case, writing uj = uidx\ A ... A dx n with 
u) G C°°(£l) strictly positive, we also notice that it is not restrictive to assume Co = 1 in fi; 
indeed, comparing ([6]) with the classical integration by parts formula in Q with no weight, 
we immediately see that 

div^X = divJT + X logo), 

where, in the right hand side, divX is the Euclidean divergence of X. One can then compare 
the integration by parts formulas in the weighted and in the classical case to obtain that 
Dxu depends on uj through the factor uj. This is also evident in the smooth case, where 
the function Xu is clearly independent of uj, but Dxu = (Xu)m. 

After this reductions to the standard Euclidean setting, we fix an even, smooth convo- 
lution kernel p in W 1 with compact support and denote by p e (x) = e~ n p(xje) the rescaled 
kernels and by u* p e the mollified functions. We shall use the so-called commutator lemma 
(see for instance [9]) which ensures 

(D x v) * p £ - X(v * p £ ) -)■ as e | 0, strongly in L} 0C (Q) (16) 

whenever v G Lj oc (£l) and X is a smooth vector field in fl. Since (Dxu) * p e m and 
\(Dxu) * p e \m converge in the duality with C c (fl) to Dxu and l-Dx^l respectively (see 
for instance |11| Theorem 2.2]), (|16p shows that the same is true for X{u * p £ ) m an d 
\X(u * p e )\m. The same proof works for vector- valued measures (convergence of total 
variations, the only thing that cannot be obtained arguing componentwise, is still covered 
by [HI Theorem 2.2]). □ 

The following result provides a characterization of in terms of difference quotients 

involving the flows generated by the vector field X (for similar results in the context of 
doubling metric spaces supporting a Poincare inequality, see [36]). 

Given a smooth vector field X in Q, C M open, we denote by the flow generated 
by X on M. By compactness, for any compact set K C £1 the flow map starting from 
K is smooth, remains in a domain £1' <g Q and is defined for every t G [— T, T], with 
T = T(K, X) > 0. Recall that the Jacobian J<3?f of the flow map x i— > $^(x) is the smooth 
function JQ^ satisfying ($^)*u; = JQ^uj, so that the change of variables formula 

J <f)UJ = J 4>o<$>fj<$>fuj 
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holds. By smoothness, there is a further constant C depending only on X (and T) such 



that J§ x satisfies 



| J$f (x)- 1| <C|t|, \J&? ( x ) - 1 -tdiv w X(x)| < Ct 2 Vx € K, t € [-T,T\. (17) 

Estimate (|17p is a simple consequence of Liouville's theorem, showing that the time deriva- 
tive of t H> log(J$f (x)) equals (div w X)($f (x)). 

Theorem 2.5. Lei C M be an open set and let u S L 1 (Q,m). Then Dxu is a signed 
measure with finite total variation in Q if and only if 

SU p | f K£f)jiH_ ^ | o < |t| < T(K,X)\ < oo. (18) 

K C H compact {J K 1^1 J 

Moreover, if Dxu is a signed measure with finite total variation in Q, it holds 

\D x u\(n) =sup(liminf f \ u ^t )~u\ u i n > m n \ ^ 
{ *-^o J n , \t\ J 

Proof. Let us prove that the existence of Dxu implies (fT8|) . For < \t\ < T(K,X) we will 
prove the more precise estimate 

r Hmz^ u < {1 + C \t\)\D xV \(a t ) with n t -.= JJ (20) 

JK ' ' r€[0,|t|] 

which also yields the inequality > in (|19p . In order to prove ()20|) we can assume with no 
loss of generality, thanks to Theorem 12. 44 that u S C°°(f2). Under this assumption, since 
the derivative of 1 1-> ii(3>^(x)) equals Ia($f (i)) we can use Fubini's theorem and (fT7]) to 
get, for i > (the case i < being similar) 

f \u(&?)-u\u< [ [ \Xu\{<f> x )udr < [ (1 + Cr) / |Xn|wdr. 
Jftr Jo jk Jo Jn t 

Estimating (1 + Cr) with (1 + Ct) we obtain (f20j) . 

Let us prove the inequality < in (|19p . By the lower semicontinuity on open sets of 
the total variation of measures under weak convergence and the inner regularity of |-Dx^|> 
it suffices to show that for all Q' (s= ft the difference quotients t _1 (n( ( 5^') — u)m weakly 
converge, in the duality with C C (Q'), to Dxu. By the upper bound (fT8|) we need only 
to check the convergence on C£°(Q') test functions. This latter convergence is a direct 
consequence of the identity (which comes from the change of variables x = & x t (y)) 

f u(<$> x )-u f y(<5> x t )J$> x t - if 

/ (pdu{x) = - udu(y) ip e C c (fi ), 

Jn> 1 Jo. 1 ~ l 

of the expansion (|17p and of the very definition of Dxu. The same argument can be used 
to show that finiteness of the supremum in f)18|) implies the existence of Dxu. 
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2.4 BV functions in metric measure spaces (X, d, m) 

Let (X, d) be a metric space, and define for / : X — > R the local Lipschitz constant (also 
called slope) by 

V/ (s) :=limsup . (21) 

y^x d(y,x) 

If we have also a reference Borel measure m in (X, d), we can define the space BV{£1, d, m) 
as follows. 

Definition 2.2. Let dclbe open and u G L 1 (Q,m). We say that u G .BV(r2,d, m) if 
there exist locally Lipschitz functions u n convergent to u in L}(Q.,m), such that 

lim sup / |Vu n | dm < +oo . 
n Jn 

Then, we define 

||Z)u||(J)) := inf < liminf / |Vu n | dm \ u n G Lip loc (fi), lim / \u n — u\ dm = > . 

I n ^°° Jn n Jn J 

In locally compact spaces, in [36] (see also |10] for more general spaces) it is proved 
that, for u G BV(S},d,m), the set function A \— > \\Du\\(A) is the restriction to open sets 
of a finite Borel measure, still denoted by | [ JOtt ] [ - Furthermore, in |36j the following inner 
regularity is proved: 

uG L 1 ^) nWi oc (ft, ci,m), sup < oo =^> u£BV(n,d,m). (22) 

n'^n 

In the next sections we shall use a fine property of sets of finite perimeters, proved within 
the metric theory in [7] (see also [6] for the Ahlfors regular case) to be sure that the set of 
"good" blow-up points has full measure with respect to ||Z) fl ls||- The basic assumptions 
on the metric measure structure needed for the validity of the result are (in local form): 

(i) a local doubling assumption, namely for all K C X compact there exist f > and 
C > such that m(B 2r (x)) < Cm(B r (x)) for all x G K and r G (0,f); 

(ii) a local Poincare inequality, namely for all K C X compact there exist f, c, A > such 
that 

/ \u — u X}T \ dm < cr / \Vu\dm (23) 
JB r (x) Jb Xt (x) 

for all u locally Lipschitz, x G K and r G (0, f), with u x , r equal to the mean value of 
u on B r (x). 

Proposition 2.6. [7] Assume that (i), (ii) above hold and let E C O be such that 1e G 
BV{n,d,m). Then 

llm|nf ml.{m(g r (x)ng) m(g r (x)\g)} > jDlj < ^ 

rio m{B r (x)) ri o h(B r (x)) 

for H-Dl^H-a.e. x G where h(B r (x)) = m(B r (x))/r. 
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In order to apply Proposition 12.61 for our blow-up analysis, we need to provide a bridge 
between the metric theory outlined above and the differential theory described in Section [2T3l 
As a matter of fact, we will see that in the setting of Section \2.3\ under mild assumptions, 
we have always an inclusion between these spaces, and a corresponding inequality between 
||J5 fl it|| and H-D^ll- First, given a function G : TM — > [0, oo] in a smooth manifold M as in 
Section [2T3l we define the associated Carnot-Caratheodory distance by 



dcc(x, y) := inf ij VG yt (jt) I T > 0, 70 = x, -y T = V j , (25) 

where the infimum runs among all absolutely continuous curves 7. Notice that since G x is 
infinite on T X M \ T>(x), automatically the minimization is restricted to horizontal curves, 
i.e. the curves 7 satisfying 7 G T>(-y) a.e. in [0, T\. We will often use, when convergence 
arguments are involved, an equivalent definition of d cc {x, y) in terms of action minimization: 

d 2 cc (x,y) :=inf ( f G Tt (7t) | 70 = x, 71 = y\ . (26) 







Notice also that we cannot expect d cc to be finite in general, hence we adopt the con- 
vention inf = +00. However, any X G r 9 (f2,D) induces, via the flow map , admissible 
curves 7 in (|25[) with speed G{^) less than 1, for initial points in Vt and \t\ sufficiently small. 
It follows immediately that for SJ'gS! and \t\ sufficiently small, depending only on Q' and 
X, it holds: 

dcc(&?(x),x)<\t\ Vx e O'. (27) 



Theorem 2.7. Let m be defined as in @ and assume that the distance d cc defined in (|25p 
is finite and induces the same topology of M . Then for any open set £1 C M we have the 
inclusion BV(Q,d cc ,m) C BV{VL,g,uj) and \\D g u\\(£l) < \\Du\\(Q). 

Proof. Take a locally Lipschitz function u in ft (with respect to d cc ) with |Vit| dm finite, 
X G T 9 (Q,V) and apply ([27]) to obtain that \u(^f(x)) — u(x)\/\t\ is uniformly bounded as 
\t\ ,|,0 on compact subsets of f2 and 

\u($? (x)) — u(x)\ ,„ , . , n 
limsup I , — < |Vit|(x) Vx G O. 

By integrating on fl' g fl we get from Theorem 12.51 that Dxu is a measure with finite total 
variation in 0, and that 



\D x u\(n) < / |Vu| dm. 
Jn 



Eventually, we apply the very definition of BV{Sl, d cc , m) to extend this inequality to all 
u G BV(£l,d cc ,m), in the form \Dxu\(£l) < \\Du\\(Q). We can now use the arbitrariness of 
X to get the conclusion. □ 
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3 Sub-Riemannian manifolds 



A frame free approach to describe sub-Riemannian structures locally generated by families 
of vector fields [37J relies on images of Euclidean vector bundles. Recall that a Euclidean 
vector bundle is a vector bundle whose fiber at a point x is equipped with a scalar product 
(•, -) x which depends smoothly on x (in particular smooth orthonormal bases locally exist). 

Definition 3.1 (Sub-Riemannian structure). A sub-Riemannian structure on ¥ is a pair 
(U, /) where U is a Euclidean vector bundle over M and / : U — > TM is a morphism of 
vector bundles (i.e., a smooth map, linear on fibers, such that f(U x ) C T X M, where U x 
denotes the fiber of U over x) such that 

Lie x (£>) = T X M, Vx G M, (28) 

where we have set 

V:={foa\aeT(U)}, (29) 
and T(U) := {a G C°°(M,U) | a(x) G U x }). 

Given a sub-Riemannian structure on M, we denote by T>(x) the vector space f(U x ) C 
T X M and we define the quadratic forms G x : T X M — > [0, oo] by 

q ( v \ _ f min {l u lx I u G U x , f(u) = v}, v£V(x) 
x[V) \+oo, v(£V(x). { ' 

Notice that in general the dimension of T>(x) need not be constant and (x,v) \— > G x (v) is 
lower semicontinuus. Let g x : T>(x) x T>(x) — > M be the unique scalar product satisfying 

9x(v,v) = G x (v) Mv G V(x). 

We shall denote by V x : T>(x) — > JJ X the linear map which associates with v the unique 
vector u G f~ l {y) having minimal norm. It can be computed intersecting f~ l (v) with the 
orthogonal to the kernel of f\\j x . 

We will often compute G in local coordinates as follows: let a\, . . . , a m be an orthonormal 
frame for (where m = rankU) in an open set O C M. Then, defining Xj = f o aj, for 
every x £ Q and v G T>(x) we have 

{m m ^ 

Y,4\v = J2 c iMx)\- (31) 
i=l i=l J 

In this case we shall also view V x as an R m -valued map. Notice that, by polarization, it 
holds 

g x (v,w) = (T x (v),T x {w)) v,we V(x). (32) 
Remark 1. The above definition includes the cases (see also Section lBTTI for more examples): 
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• U is a subbundle of TM and / is the inclusion. In this case the distribution T> has 
constant rank, i.e., dim.T>{x) = dimU^ = rankU. When U = TM and / is the 
identity, we recover the definition of Riemannian manifold. 

• U is the trivial bundle of rank m on M, i.e., U is isomorphic to M x M m and T> 
is globally generated by m vector fields / o e±, . . . , / o e^, where ej(x) = (x, Ej) and 
e±, . . . , e m is the canonical basis of M. m ; in particular, we recover the case when M = 
f2 C R n and we take m vector fields satisfying the Hbrmander condition. 

The finiteness of the Carnot-Caratheodory distance d(-, •) induced by G as in (|25p (note 
that we will drop from now on the cc in (|25p and (|26p ) is guaranteed by the Lie bracket 
generating assumption on D (see [5]), as well as the fact that d induces the topology of M 
as differentiable manifold. The metric space (M, d) is called a Carnot-Caratheodory space. 

Since T> is Lie bracket generating, at every point x G M there exists k x G N such that 
the flag at x associated with T> stabilizes (with step k x ), that is, 

{0} C V 1 (x) C V 2 (x) C • • • C V k *(x) = T X M, (33) 

where P 1 (x) = T>(x) and T> t+l {x) = (T) % + [D, T> l ])(x). The minimum integer k x such 
that ()33p holds is called degree of non-holonomy at x. With the flag (|33h we associate two 
nondecreasing sequences of integers defined as follows. The growth vector is the sequence 
(ni(x), . . . ,nk x (x)), where rii(x) = dimP l (x). Notice that for every x G M, rik x {x) = n. To 
define the second sequence, let v±, . . . , v n G T X M be a basis of T X M linearly adapted to the 
flag ([33]) . The vector of weights is the sequence (w± (x) , . . . ,w n (x)) defined by Wj(x) = s 
if Vj G T> s (x) \ T> s ~ 1 (x). Notice that this definition does not depend on the choice of the 
adapted basis, and that 1 = W\(x) < ■ ■ ■ <w n {x) = k x . 

We say that a point x G M is regular if the growth vector is constant in a neighborhood 
of x, otherwise we say that x is singular. If every point is regular, we say that the sub- 
Riemannian manifold is equiregular. 

3.1 Examples 

In this section we mention some examples of sub- Riemannian manifolds. A first fundamental 
class of examples is provided by Car not groups. 

Example 3.1 (Carnot groups). Let us consider a connected, simply connected and nilpotent 
Lie group G, whose Lie algebra q admits a step s stratification 

q = V 1 ®V 2 ® ■■■ ®v s , 

namely [Vi, Vj] = Vj + \ for every j = 1, . . . , s — 1 and [V±, V s ] = {0}. Every layer Vj of g 
defines at each point x G G the following fiber of degree j at x 

W x = {Y(x) G T X G | Y G Vj} . 

All fibers of degree j are collected into a subbundle H 3 of TG for every j = 1, . . . , s. Fix 
a scalar product (•, -) e on H\. Then this canonically defines a scalar product (•, -) x on H\ 
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by left invariance. In this way we endow the vector bundle H with a Euclidean structure. 
In the language of Definition 13, 1\ the inclusion i : H 1 — > TG defines a left invariant sub- 
Riemannian structure on G. According to (|29|) . the corresponding module T> is precisely 
made by all smooth sections of H , the so-called horizontal vector fields and T>{x) = H% for 
every x E G. The validity of (|28p is ensured by the assumption that is stratified. Note 
that different choices of scalar product on H\ define Lipschitz equivalent sub-Riemannian 
structures on G. 

The group G equipped with a left invariant sub-Riemannian structure is called Carnot 
group. Let m be the dimension of Vi, n the dimension of G, and let (f/i, . . . , y n ) be a system 
of graded coordinates on G. An equivalent way to define a left invariant sub-Riemannian 
structure on G is the following. Fix a basis X\ , . . . , X m of V\ with the following form 

n 

Xj(y) = dj + ^2 a ^(y) di for ever y 3 = 1 ' ■ ■ ■ » m ' 

i=m+l 

where Qjji IS El homogeneous polynomial such that aji(5 r y) = r uli 1 aji(y) for every y £ G 
and r > 0, where <5 r y = Y^=i rUlj ^j and the degree Uj is defined by the condition ej G V^^ 
for every j = l,...,n. Under these coordinates, we take the Euclidean vector bundle 
U = GxK m and 

in 

f : G x R m ->■ TG, /(y,0 = (y,E^' X ^y))- 

Example 3.2 (Heisenberg group). The Heisenberg group is a special instance of a step 2 
Carnot group. It can be represented in the language of Definition 13.11 as ]R 3 equipped with 
the left invariant vector fields 

x i=di~ y<9 3 , X 2 = d 2 + y<9 3 , 

with respect to some polynomial group operation. We have V\ = span{Ai,A2} C f) and 
V<2 = span{A3} C f), where f) is the 3-dimensional Lie algebra of left invariant vector fields 
and A3 = ^3. Following the previous general case of Carnot groups, we set U = M 3 x M 2 , 
equipped by the morphism / defined by 

f(x,(l,0)) = X 1 (x), /(x,(0,l)) = A 2 (x), 

for every x G M 3 . In this case, we have T> = {b\X\ + b 2 X 2 : b\ , b 2 : M 3 — > R smooth}. 

Example 3.3. Let us consider a connected and nilpotent Lie group G with Lie algebra 53 
equipped with a linear subspace q± = span{Ai, . . . , X m } that satisfies 

= 01 ©02 and [0i,0i] =02- 

We notice that here need not be nilpotent and G is not necessarily simply connected, 
hence the exponential mapping in general is not invertible. This cannot occur for Carnot 
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groups, where the exponential mapping is always bianalytic. The main point is that the 
groups G in general need not have dilations. Here the Euclidean vector bundle defining 
their sub-Riemannian structure is given by 



/ : G x R m -»■ TG, f(y,0 = {y,Y,Zj X j(y))- 

i=i 

The foremost example of these groups is the rototranslation group M 2 x S , m = 2, equipped 
with vector the fields 



X\ = cos6d x + sin 9d y and X 2 = dg. 

In the next examples, we refer the reader to Sections 13.31 13.41 for notions of nilpotent 
approximations and privileged coordinates. 

Example 3.4 (Rank- varying sub-Riemannian structure - Grushin plane). Consider the 



and /((xi,x 2 ),(l,0)) = X h 



sub-Riemannian structure on R 2 defined by U = R 2 x 
f((x 1 ,x 2 ), (0, 1)) = X 2 , where 

Xi = di, X 2 = x\d 2 . 

Then V(x) = spanjX^iz;), X 2 (x)} and n 1 (x 1 ,x 2 ) = 2 if x\ 7^ 0, whereas rii(0, x 2 ) = 1. The 
growth vector at points in the vertical axis E = {(xi,x 2 ) \ x\ = 0} is equal to (1,2). On 
the other hand, at points of M 2 \ E, the growth vector is equal to (2). In other words, E is 
the set of singular points and the varying dimension is n\{x). Given v = v\d\ + v 2 d 2 , the 
sub-Riemannian metric in this case is 



G x (v) 



x i 

xi = 0,v 2 



00, 



As a consequence, the scalar product is 



g x (v,w) = v\wi + 



v 2 w 2 







x\ = 0, v 2 7^ 0. 



v,w*EV(x), v = v\di + v 2 d 2 , 



w 



wid\ + w 2 d 2 . 



At (0, 0) (and at any point in the vertical axis), the nilpotent approximation (and thus the 
metric tangent cone) is the sub-Riemannian structure itself. Indeed, since the degree of 
non-holonomy at (0,0) is 2 and since coordinates (xi,x 2 ) are linearly adapted to the flag 
of V at (0,0), they are also privileged. The weights at (0,0) are wi(0,0) = 1,^2(0,0) = 2. 
Hence both X%, X 2 are homogeneous of non-holonomic order —1 at (0, 0). Geodesies can be 
computed explicitly and the Carnot-Caratheodory distance is homogeneous with respect to 
the dilation ^(^1^2) = (Axi,A 2 X2). 
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Example 3.5 (Singular point at which the metric tangent cone is a Carnot group). Consider 
the sub-Riemannian structure on R 3 given by U = K 3 x K 3 and /(x, (1,0,0)) = X\, 
f(x, (0, 1, 0)) = X 2 , f(x, (0, 0, 1)) = X 3 , where 

x 1 = a 1 -^-a 3 , x 2 = d 2 + ^d 3 , x 3 = x 2 3 d 3 . 

Set V(x) = span{Xi(x), X 2 (x), X^(x)}. Then ni(xi, x 2 , x%) = 3 if X3 7^ 0, whereas 
ni(xi, x 2 , 0) = 2. The growth vector at points in the plane £ = {(xi,X2,X3) | X3 = 0} 
is equal to (2,3). On the other hand, at points of M 3 \ S, the growth vector is equal to 
(3). In other words, X is the set of singular points and the varying dimension is n\. The 
sub-Riemannian metric is 



G x (v) 



vf + X 3 = 0, 2u 3 + X 2 ^l - X^ = 

00, x 3 = 0, 2u 3 + x 2 vi - x\v 2 j=- 0. 



Coordinates (xi, £2, £3) are linearly adapted to the flag of 2? at (0,0,0) and the degree of 
non-holonomy of the structure at (0,0,0) is 2. Hence {xi,x 2 ,x%) are privileged at (0,0,0). 
A simple computation shows that ordo^i = ordo-^2 = — 1 whereas ordoX3 = 1. Therefore, 
the truncated vector fields are X\ = X\ , X 2 = X 2 and X3 = and the metric tangent cone 
at (0,0,0) is isometric to the Heisenberg group (see Example 13. 2p . 

Example 3.6 (Generalized Grushin plane). Let us consider a generalization of Example 13.41 
where we replace X 2 with 

X%(x)=x 2 *d 2 , 
with a > 1. The sub-Ricmannian metric becomes 



x l 

vf, xi = 0, v 2 = 

+00, xi = 0, v 2 i=- 0, 

from which we deduce that the map V x '■ T>{x) — > M. 2 is given by 



V x (v) 



(t>i,0) xi = 0, u 2 = 0. 



The growth vector is (1, 2) at points in the vertical axis and it is (2) outside the vertical axis. 
Set X(x\,x 2 ) = x\d 2 . Then X(x) £ T>(x) for every point x but x 1— )• Ga;(X(x)) explodes at 
points in the vertical axis. Indeed, 



V x {X{x)) 



0,^t), xx t^O 
(0,0) xi =0, 



IS 



whence 

10 X! = 0. 



G x {X{x)) = \V x {X{x))Y 



Notice however that x h- > ^(^(x)) is measurable. 

The sub-Riemannian structures of examples l3.4l[33ll3.6l are also called almost-Riemannian, 
see OH]. 

Example 3.7 (corank-1 or contact distributions [26]). Let M be a smooth manifold and 
/3 be a completely non-integrable one-form on M. Set U = ker/3. Then U is a vector 
bundle in M of rank dimM— 1. Choosing any Euclidean structure on U, we can define the 
sub-Riemannian structure (U, i) on M where i is the inclusion. The growth vector of the 
distribution is constantly equal to (dim M — 1, dim M) and the structure is equiregular. This 
class of sub-Riemannian manifolds satisfies at each point the assumptions of our blow-up 
theorem below, see Section UJ 

3.2 BV functions on sub-Riemannian manifolds 

In this section we provide characterizations for BV functions in sub-Riemannian manifolds 
and prove the Riesz theorem. First of all, we notice that in a sub-Riemannian manifold one 
can locally fix an orthonormal frame 

X\ = f o at, X 2 = f o ct 2 , • • • , X m = f o a m (34) 



where a±, . . . ,a m is a local orthonormal frame of U. The frame (JMJ) defines the vector 
measure 

X.u := {D Xl u, . . . ,D Xm u). (35) 

Theorem 3.1. Let Q, C M be an open set and let u G L l {Sl). Then, the following three 
conditions are equivalent: 

(i) sup{\D x u\(n) I X = foa, ae r(U| n ), M < 1} < oo; 

(ii) u G BV(Q,d,m); 

(Hi) u G BV(fl, g, uj). 

Furthermore, if one of the previous conditions holds, then we have 

\\D g u\\(Q) = \\Du\\(n) = sup{\D x u\(n) \ X = foa, a £ T(\J\ U ), \a\ < 1} . 

If Q has an orthonormal frame (|34j) . then \\D g u\\(£l) = |Xu|(f]), where Xu is the vector 
measure ()35[) defined on Q. 
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Proof. For every open set A C 0, we define the set function 

s(A) = sup {\D x u\(A) \ X = f o a, aeT(V\ A ), \a\ < 1}. 

The simple inequality s(Cl) < \\D g u\\(Cl) is a consequence of the fact that a larger class of 
vector fields is considered in the definition of ||D 9 u||, hence the implication from (iii) to (i) 
follows. From Theorem 12.71 we get ||D 9 u||(r2) < ||Z)u||(J)), hence the implication from (ii) 
to (iii) follows. Next, we prove the implication from (i) to (ii), that follows by establishing 
the inequality 

\\Du\\{n) < s(Q). 

By a partition of unity, we can assume with no loss of generality that in the vector fields 
Xi = f o o~i are globally given. We assume first that u € C 1 (ri). In this case we prove first 
the inequality (where the left hand side should be understood as the slope (|2ip w.r.t. d), 
in local coordinates 



\Vu\ 2 (x) < ^2(X iU (x)) 2 . 



(36) 



i=l 



In order to prove this inequality, if c £ L 2 ([0, l];R m ), 7 = X^ c «^«(7)> To = x an d 71 = y, 
we have 



\u(x) - u(y)\ 



d lt u(j t 



„\ m m 

dt \ = \ y2ci(t)Xiu(-yt)dt\ < \\c\\ 2 sup V(X, 
J o i= i te[o,i] \ i=1 



Minimizing with respect to c gives 
\u(x) - u(y)\ 



d{x,y) 



< 



sup 



J2(X iU (z)) 2 I d(x,z)<2d(x,y)\. 



Then, taking the limit as y — > x provides (|36p . Now, considering the vector- valued measure 
Xu in (|35|) . whose total variation |Xu| is equal to ^^.^XiU^m, we may write 

/ |Vn| dm < \Xu\(Q). 

We can now invoke the definition of BV(£l, d, m) and Theorem 12.41 to obtain the inequality 

||Du||(n) < |X«|(n). (37) 

The estimate |Xu|(f2) < s(Q) immediately follows observing that for each (p £ C,J(r2,]R m ) 
with \ip\ < 1 there holds 

„ m 

/ u div w ( y^y?j Xj) u < s(fi) . 
Jo i=1 

Collecting all previous inequalities, we achieve 

s(n) < \\D g u\\(n) < \\Du\\(n) < \Xu\(n) < s (n), 
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that establishes all the claimed equalities. □ 
From now on, in view of Theorem 13 .1\ the measures and \\Du\\ will be identified, 

and we will use their local representation as total variation of the vector- valued measure Xu 
in ([33]) . We also notice that due to Proposition 12.21 and the Poincare inequality with respect 
to vector fields in R n , see for instance [30] and [32], a local Poincare inequality also holds in 
our framework. This implies that we can apply Proposition 12.61 to obtain the inequalities 

r|o m{B r {x)) ri0 h(B r (x)) 

(recall that h(B r (x)) = m(B r {x))/r) for ||£? g l^||-a.e. i£!!, whenever 1 E £ BV(fl, g,uj). 

Definition 3.2 (Dual normal and reduced boundary). Write, in polar decomposition, 
Xl# = z/gll-Dgl^H, where u% : f2 — > W 71 is a Borel vector field with unit norm. We 
call v E dual normal to E. 

We denote by J^E the reduced boundary of E, i.e. the set of all points x in the support of 
1 1 D g 1 E \ I satisfying ([38]) and 

lift -i iL , ^ / |^(y)-^(x)| 2 ci|| J D 3 l E |Ky) = 0. (39) 

It is simple to check that, while the dual normal v* E depends on the choice of the 
orthonormal frame, the reduced boundary J~gE does not. 

We notice that (I38p and the relative isoperimetric inequality give 



. t \\Dgl E \\(B r {x)) \\D g l E \\(B r (x)) ^ 

< lim inf — \ < hm sup — " < oo (40) 

rio h(B r (x)) ri0 h{B r (x)) 

Then, the doubling property of h implies the asymptotic doubling property: 

\\D g l E \\{B 2r {x)) 
hmsup n — - — rrj— . < oo tor \\D g i-E -a.e. x G \l. 
r\.0 \\Dgl E \\{B r (x)) 

We shall use the following proposition, a direct consequence of the Lebesgue continuity 
theorem in all metric measure spaces with an asymptotically doubling measure: here our 
measure is 

Proposition 3.2. // E has locally finite perimeter in £1, then \\D g l E \\-a.e. point of Q 
belongs to J~gE. 

Now, we are in the position to establish Riesz theorem in sub-Riemannian manifolds, 
compare with Remark [3j As a byproduct, applying Riesz theorem to a characteristic func- 
tion 1 E , we can identify a geometric normal u E , image of the dual normal under the mor- 
phism /. 
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Theorem 3.3 (Riesz theorem in sub-Riemannian manifolds). Let u G BV(fl, g, uj) . There 
exists a Borel vector field u u satisfying G{u u ) = 1 || Z? 5 ti|| -a. e. in Q. and 

D x u=g(X,v u )\\D g u\\ VlGr»(n,D). (41) 

If E is a set of finite perimeter and u = 1e, ve '■= V \ E is given in a local frame Xi = f o o~i 
by f(%2i v *e i a i) an d ^ will be called geometric normal. 

Proof. By a partition of unity, we can assume that in f2 an image of an orthonormal frame 
X\ = /oo"i, . . . ,X m = f oa m is globally given and, taking into account Proposition 12.21 we 
can assume with no loss of generality that f2 C M n and that ui = ujdx\ A ... A dx n . Let Xti 
be as in (|35p and write, in polar decomposition, Xu = io|Xu| = for some Borel 

w : £1 — > M. m with =1. If X = C{X{ with q smooth, obviously 

m m 

D x u = ^2ciDx l u = ^2c i w i \\D g u\\. (42) 

i=l i=l 

Assume now that X € T g (£l,T>). By a measurable selection theorem we can write X = 
Yli c iXi with Ci Borel and ^ c? < 1. If we define, as in the proof of Theorem 12.41 X € = 
YliLi c i * PeXi, the commutator theorem and the distributional identity div w Y = divY + 
yioga) give (notice indeed that div^QA",) makes sense only as a distribution) 

v t = div w (q * p e Xi) - [diVu(ciXi)] * p e -)• strongly in Ll oc (Q) Vi = 1, . . . , to. (43) 

Hence, adding with respect to i gives div^X 6 — > div^X strongly in Lj oc (£l). Recall that the 
distribution div(ciXj) satisfies 

f 1>(y)ddiv(ciXi)(y) = - [ Ci(y)X t i>(y)dy G ^(fl), 
Jn Jn 

whence 

div^aXi) * p e (x) = / p e {x - y)dd\\(ciXi){y) + / p e (x - y)a(y)X "Jog Q dy. 
Jn Jn 

Thus, a direct computation shows that 




where X^ are smooth functions such that Xi = Y^j=i-^-ij^j- The nrs ^ summand is locally 
uniformly bounded since \a\ < 1 and Xy is smooth. As for the second term, changing 
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variable we obtain 



/n 
Ci(y)^2djPe(x - y)(X ij (x) - Xijiyj) 
i=i 

n 1 „ 

< WcidjpW^- J \Xij(x - ey) - X i:j {x)\dy . 



3=1 

It follows that u e is locally uniformly bounded, and the convergence in (1431) holds in the 
weak* sense in L^ c as well. Therefore, up to subsequences, div^X 6 — > div^X weakly* in 
Lfo c (p,). Since u G L l {0), we can pass to the limit into (|42|) with X = X e to obtain that 
(|42p holds for any smooth vector field X = c^X with Cj just bounded Borel. 

Now, let us prove (f4Tj) with v u = f(w) = Yli w i^i- To this aim, we notice that in the 
representation X = Y^i CiXi we can always assume that c(y) is orthogonal to the kernel of 
f\u . Recalling that 

m mm 



aA = 9y(%2 aiXi(y), ^ kXi(y)) 



for all a orthogonal to the kernel of f\u y and all b € M m , we apply the previous equality 
with a = c and b = w, to obtain (|41|) . Notice that it has been essential the estabilishment 
of (I42p with non-smooth c's: even if the initial c's were smooth, their pointwise projection 
on the orthogonal to the kernel of / might be not smooth (see Example I3.6|) . 

By construction, G{y u ) < 1, because |io| = 1; the converse inequality can be proved 
noticing that on any Borel set A it holds \D x u\(A) < f A ^G{v u ) d\\D g u\\, for all X G 
T g (£l,T>). Choosing A open and maximizing with respect to X gives 



|1V||GA)< I VGWjd\\D g u\\. 

J A 



Since A is arbitrary, we have y G(u u ) > 1 ||D g u||-a.e. in Q. □ 

Remark 2. A byproduct of the previous proof (just take w = v 3 ^, in the previous proof, 
and notice that we proved that G(J2 i WiXi) = 1) is the fact that the dual normal v* E is 
orthogonal ||Z) fl l^;||-a.e. to the kernel of /. 

Remark 3. It is rather natural to ask whether Theorem 13.31 holds in the general framework 
of Section [2.31 where G is only Borel and we consider general smooth sections of T>. More 
precisely, in this setting, for u G BV(£l, g, u), it would be interesting to find a positive finite 
measure a in Q and a Borel vector field v u in Q with G(v u ) = 1 ||J5 fl it||-a.e., satisfying: 

D x u = g{X,u u )a VX G r»(fi,X>). (44) 

A good candidate for the measure a would be the supremum of \Dxu\, in the lattice of 
measures, as X varies in T g (Q,T>). 
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3.3 Privileged coordinates and nilpotent approximation 

In this section we recall the notion of privileged coordinates and of nilpotent approximation 
of a sub-Riemannian manifold at a point. 

Let (U, /) be a sub-Riemannian structure on M and fixp£ M. Let Q be a neighborhood 
of p and let . . . , a m be a local orthonormal frame for the U|q, where m = rankU. Define 
Xj = foaj. 

Given a function ip G C°°(M), for i G {l,...,m} we call X^ip a first non-holonomic 
derivative of ip. Similarly, if i, j G {1, ...,m}, XiXjip is a non-holonomic derivative of 
order 2. With this terminology, we say that ip has non-holonomic order at p greater than 
s if all non-holonomic derivatives of ip of order a < s — 1 vanish at p. If moreover there 
exists a non-holonomic derivative of order s of ip which does not vanish at p we say that ip 
has non-holonomic order s at p. 

By duality, given a differential operator Q, we say that Q has non-holonomic order > s 
at p if Qtp has order > s + n at p whenever ip G C°°(M) has order > 77 at p. Clearly, the 
non-holonomic order (of a function or of a differential operator) is an intrinsic object, i.e., 
it does not depend on the chosen vector fields X%, . . . ,X m . 

Definition 3.3 (Privileged coordinates). Let (p = (ip\, . . . ,f n ) '■ ^ - > ^ n be a coordinate 
system centered at p, i.e., (p is a smooth diffeomorphism and <p(p) = 0. We say that (p is a 
system of privileged coordinates if 

• the canonical basis (d Zl , . . . , d Zn ) of TqW 1 is linearly adapted to the flag associated 
with 93* D at 0; 

• for every i = 1, . . . , n the non-holonomic order of the i-th coordinate function z 1— >• z% 
at is equal to Wi{p). 

Existence of privileged coordinates at points of sub-Riemannian manifolds have been 
proved in a constructive way in several works |15 t l27 1 l29 | l39j. Moreover, if the non-holonomy 
degree k p is 2 at a point p (see Example 13. 7p . each coordinate system satisfying the first 
property in Definition 13.31 directly satisfies the second one. 

Let ip : Q — > M. n be a system of privileged coordinates at p. We consider the sub- 
Riemannian structure (U|^, <p* f) on W 1 . Clearly, the vector fields (p*Xi, . . . , ip x X m are 
global generators for ip*T>. Using Proposition 12.2^ the order of a function tp G C°°(£l) at p 
coincides with the order of ip o ip" 1 G C°°(]R n ) at 0. 

Privileged coordinates allow to compute non-holonomic orders (both of functions and 
of differential operators) using the following facts. 

(i) A monomial function h G C°°(lR n ), h(z) = z^z^ 2 ■ ■ ■ z% n has order w\{p)ai + • • • + 

VJn(p)®n at 0. 

(ii) Given i G {1, . . . , n}, a vector field F(z) = z^z^ 2 • • • z® n d Zi has order w\{p)a\ + • • • + 
Wn(p)a n ~ Wiip) at 0. 
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Thanks to (i), the order at of a function h G C 00 (]R ri ), denoted with ordo(/i) is the 
smallest number w\{p)a\ + • • • + w n (p)a n , such that a monomial z® 1 z% 2 ■ • • z" n appears 
with a nonzero coefficient in the Taylor expansion of h at 0. Using (ii), we have a notion of 
homogeneity of vector fields. Namely, a vector field F on M. n is homogeneous of order s if 

n 
i=l 

where 

ord (/») - Wi(p) = s, Vz = l,...,n. 
By definition, the order of (p*Xi at is greater than —1. Hence, we have an expansion 

(s) 

where K is the homogeneous component (p*Xi of order s. Define m vector fields on W l 
by 

Xi = Yt X) - (45) 

Denote by V the distribution on W 1 generated pointwise by X\, . . . ,X m and define, in 
analogy with (|3T]l 

min { E cf | v = Yl CiXi(x) > , v £ V(x) 

U=i i=\ ) _ (46) 

+oo, v <£ V(x), 

and g x the corresponding scalar product on T>(x). 

Remark 4. Take U = W 1 x R m and / : U TR n defined by f(z,v) = TT=i v A z )- 
Then one can define T> and G x as the one induced by the sub-Riemannian structure (U, /) 
on R n . The fact that Lie^P = W 1 for every z € M n follows by the Lie bracket generating 
condition on ip^V. 

Denote by d the Carnot-Caratheodory distance on R n associated with the sub-Riemannian 
structure (U, /), and denote by B r the set {y E W l \ d(y, 0) < r}. Given A > 0, define the 
dilation 5 X : M. n -»• M n by 

5 x ( Zl , ...,z n ) = (X w ^ Zl , X w ^z n ) (47) 

The sub-Riemannian structure (U, /) on IR n defined in Remark [3] is called a nilpotent 
approximation of (U, /) at p. Let us recall some properties of nilpotent approximations 
that will be useful in the sequel. 

Proposition 3.4. Let (U, /) be a nilpotent approximation of (U, /) at p. Then: 
(i) the growth vector of (p*T> at coincides with the growth vector of T> at 0; 
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(ii) any vector field V G Lie{Xi, . . . , X m } is complete and Lie{Xi, . . . , X m } is nilpotent; 

(Hi) the distance d is homogeneous with respect to 5\, i.e., d(5\z,5\z') = Xd(z,z'), for 
every X>0,z,z'g R n ; 

(iv) given r > and a smooth vector field X on Q such that ord p X > — 1, the vector field 
Y r on W 1 defined by 

Y r = r{5 1/ Mv*X-X), 

where X is the homogeneous component of ip^X of order —1 atO, satisfies the following 
property: Y r and its divergence converge uniformly to zero on compact sets of W 1 as 
r tends to zero. 

Proof. It is easy to see that if X has order > a and Y has order > j3 at then [X, Y] has 
order > a + (5 at 0. If X is homogeneous of order a at and Y is homogeneous of order f3 at 
then [X, Y] is homogeneous of order a + ft or it is zero. Let Xj = [X- Lk [• • • [Xi 2 , X- Ll \ •••]], 
where I = . . . , i^) G {1, . . . m} k . Denote by Xj the Lie bracket [Xi k [• • • [Xi 2 , X^] ■■■]]■ 
Since Xj is homogeneous of order —k (or it is zero), (p*Xj — Xj has order > — k at 0. 
Therefore 9?*X/(0) - X/(0) G ip*V k ~ l (0). As a consequence, dimP fc (0) = dimip*V k (0) for 
every k, which gives the first property. 
By homogeneity, for every i = 1, . . . , m, 

n d 

j=l 3 

with fij satisfying 

f ij (5 x z)=\ w i- 1 f ij (z). 

This implies that fij is a homogeneous polynomial of non-holonomic degree Wj — 1, whence 
it depends only on coordinates z^ with k such that Wj-(p) < Wj(p). Let j < ni (where 
ni = dimD(O)). Then, since fij is constant, the solution of Zj = fij(z) is a linear function 
of t. Take now j G {n! + 1, ...,712} then, since /^(z) only depends on Z\, . . . , z ni , the 
solution of = fij(z) is a quadratic function of t. Iterating this process we obtain that the 
flow of Xi is defined for every t, that is, Xi is complete. 

Since fij is a polynomial of degree Wj{p) — 1 < w n (p) — 1 = k p — 1, every Lie bracket 
between the Xj of length greater than /c p vanishes identically. Therefore, the Lie algebra 
generated by X\, . . . , X m is nilpotent and, for every V G Lie{Xi, . . . , X m }, we have 

n 

v = 1 £v j (z)d Zj , 

where Vj is a polynomial of non-holonomic degree < Wj(p) — 1. Using the above argument, 
one infers that the flow of V is defined for every t. 

The homogeneity of d is a consequence of the fact that, under the action of 5\, the 
length of a curve (calculated with G) is multiplied by A, which in turn follows by Xj being 
homogeneous of order —1. 
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If X is a smooth vector field on £1 having order at p greater than —1 then ip*X = X + R, 
with X homogeneous of order —1 at and R having order > at 0. Homogeneity of order 
— 1 at means that X satisfies 

(6 1/r )*X(z) = r-'Xiz), 

whence 

r{5 1/r )S(z) = X{z). 

Let R = Ya=i c i{ z )®zi- Since R has positive order, there exist po > and Co > such that 
\(H(z)\ < C (|zi| 1/Wl + • • • + \z n \ l / w -) w \ Vz e B po . (48) 

Let K C M" be any compact set and let Co = ((p" 1 )*^ (see Proposition I2.2|) . Denote by 
Co the density of Co with respect to the Lebesgue measure, i.e., Co = Codz. Thanks to the 
identity div^Y = divY + Ylogw, it suffices to show that the Euclidean divergence of Y r , 
i.e., the divergence with respect to the Lebesgue measure, converges to zero on compact 
sets. We have 



Y r (z) = r[(6 1/r %R](z)=rY / r~ m c l (5 r z)d Zi , 
div Y r (z) = 



i=i 



where hi(r, z) = Ci(S r z). Hence, to prove the required convergences, it suffices to show that, 
for every i = 1, . . . ,n, 



lim sup sup I Ci (6 r z) I < oo 

r|0 r w i K 



lim sup sup 

rlO r w * K 



dhi 



dzi 



(r, z) < oo. 



(49) 
(50) 



Assume r < po/di&mK. Then b r z E B po , whence (|48|) implies 

\ci(S r z)\ < CV i (N 1M + ■■■ + l^l 1 ^")^ < C r Wi max(|zi| 1 / , ° 1 + • • • + \z n \ 1 ^) 1 

and is proved. As for (foTIj) . we have 

dhi 



whence 



(r,,)=^( M , 



1 



lim sup sup 

rlO r w * K 



dhi 


(r, z) 


< 


dci 


dzi 




dzi 



(0), 



since c,- is smooth. 



□ 
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3.4 Nilpotent approximation and metric tangent cones 

The following theorem provides an estimate between the sub-Riemannian distance d and 
the distance d associated with the nilpotent approximation of (U, /) at p. It has been 
proved for equiregular sub-Riemannian manifolds in [34, Proposition 4.4] (see also [15] for 
the general case). Our proof is inspired by the arguments in [IJ Lemma 8.46, Theorem 8.49]. 
Just for notational simplicity, we omit the diffeomorphism (p and rename the vector fields 
ip*Xi by Xi. 

Theorem 3.5. Let d and d be the Carnot-Caratheodory distances associated with the family 
of vector fields X\, . . . , X m and Xi, . . . , X m , respectively. Let K r be the closure of B r . Then, 
the following estimate holds 

lim- sup \d(x,y) - d(x,y)\ = Vi? > 0. (51) 

Proof. If x, y € Kr^, then we write x = 5 e x, y = 5 t y, with x, y E Kr, where Kr = Si^Kr^. 
Using homogeneity of d (and renaming x,y), (f5T|) can be restated as 

d(5 e x, 5 e y) 



lim sup 

40 x , y<=K R 



d(x,y) 



\/R > 0. 



Set 

d t (x,y) = d -^^, 

and 

x i = t{h/e)* x u i = l,...,m. 

Using the last statement in Proposition 13.41 Xf converges to Xi uniformly on compact 
sets. By construction, d e is the Carnot-Caratheodory distance associated with X\, . . . , X^. 
Recall that in [1] Formula 8.26] it has been shown the existence of a constant C depending 
only the blowup point (i.e., the origin) and on the compact set Kr such that if e is small 
enough 

d t {x,y) <C\x-y\ l l k », Vx,yeK R , (52) 

where ko is the non-holonomy degree of the sub-Riemannian structure at 0. 

Since (j52l) provides equicontinuity, it suffices to show that d e — > d pointwise on Kr x Kr. 
We prove first the lim sup inequality. Set U = L 2 ([0, l],R m ) and choose c € U such that 
(recall the formulation (j26h in terms of action minimization) d(x,y) = \\c\\2 and 7(1) = y, 
where 7(0) = x and 7 = Yli ^^)- Then, if y e = 7 e (l), where 

m 

i=l 

standard ODE theory and the uniform convergence of Xf to X{ on compact sets give y e — > y. 
On the other hand, the very definition of d e gives d e (x,y e ) < 1 1 c 1 1 2 - By (I52j) we obtain that 
limsup e d e {x, y) < d(x,y). 
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In order to prove the lim inf inequality fix a sequence (e^) \, on which the lim inf e d € (x, y), 
that we already know to be finite, is achieved. Choosing c h G U such that 

rn 

j h = ^2^X t l h ( lh ), 7 h(0) = x, 7ft(l)=y, \\c h \\ 2 = d €h (x,y), 

i=l 

we can assume with no loss of generality that c h weakly converge in IA to some c. Again, 
standard ODE theory and the uniform convergence of X\ to Xi on compact sets give 

m 

j = ^2 ^iW) 7(0) = x, 7(1) = y. 
i=l 

Hence, d(x,y) < ||c||2- Since ||c||2 < lim inf ^ \\c h \\2 we obtain the liminf inequality. □ 

Remark 5. Notice that by the Ball-Box theorem (see [37]), there exists a constant L > 
such that 

B e/L (Q) C B e (0) C B Le (0) (53) 
for all e > sufficiently small. Hence, (|51jl is equivalent to 

lim- sup \d(x , y) — d(x , y)\ = MR > 0. 

The main consequence of (|5ip is that (M n ,d) is a metric tangent cone in Gromov's 
sense (see |28| ) to (M,d) at p, the quasi-isometry being the identity map (in privileged 
coordinates centered at p). Note that by very definition, a metric tangent cone carries also 
a homogeneous structure, relying on a 1-parameter group of dilations. 

Under an additional assumption, the nilpotent approximation (and thus a metric tangent 
cone) is a Carnot group. To see this, let Q be the group of diffeomorphisms of R n generated 
by the se1@ 

Uf; 1 o o • • • o ^\ ti E R,ij e {1, . . . ,m}, k G n} , 

where, obviously, we take the composition as the group operation. Thanks to the Baker- 
Campbell-Hausdorff formula, since Lie{Xi, . . . , X m } is nilpotent, for every & £ G there 
exists V G Lie{Xi, . . . , X m } such that $ = Definci 

a p = {$ G a I $(0) = 0}. (54) 

Proposition 3.6. If Q p = {Id^} then there exists a group operation * on W 1 such that Xi 
are left invariant vector fields. 

2 Recall that <f>Y denotes the flow generated by a vector field Y. 

3 We emphasize the dependence of Q p on the point p at which the nilpotent approximation is considered. 
(Recall that (f(p) = and the vector fields X\, . . . , X m actually depend on p.) 
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Proof. Define the map ^ : Q p — > W 1 by = $(0). Since the Xj are bracket generating, 

is surjective and, by assumption, \& is injective. Thus, for every x G M n there exists a 
unique $ G £ p such that $(0) = x. Taking V G Lie{X 1; . . . ,X m } such that <J> = ^f, we 
have ^^(O) = x. Notice that V may not be unique. Define the operation * : W 1 x M n — » M n 



z*y := S^oS^O). 

where V G Lie{Xi, . . . , X m }, respectively W G Lie{Xi, . . . , X m }, is a vector field such that 
^(O) = x, respectively, $1 (0) = y. Let us verify that x*y is well-defined, i.e., it does not 
depend on the choice of V and W. Let V, W be such that $r( ) = x > ®T (°) = 2/- Then > 
using ($f J" 1 = 



Thus, our assumption implies that $ x o 
z G M n . Then 



d>r w 'oci>r(o) = o. 

W „ _ TJ _ n; that iSj $W( z j = $f ( z ) for every 



$f o (0) = $r (ac) = $r' (x) =<o (0) 



It is easily seen that (R n ,*) is a Lie group, where the inverse of x = (0) is given by 
x^ 1 = &i V (0). Let : R n — > W n be the left translation, i.e., l x y = x*y. Then, by 
definition of push-forward, 



d 
~d~t 



iUiit))) 



t=o 



where j(t) = & (l x 1 y). For every t, since \I/ is bijective, there exists Z(t) G Lie{Xi , . . . , X m } 



such that $f W (0) = 7(f). Let V and W be such that ^(O) = x and $f(0) = y. We have 

*fW( ) = 7 (t) = ®hi- l y) = ^i^+y) = o $f o $ r ^(0). 



Hence, since ^ is injective, ^ 



Z(t) 



o<l>j[ 1/ o<I) 1 l/ as diffeomorphisms. Thus 



d 



t=o 



_d 
dt 
d_ 

dt 

d 
dt 



t=Q 



t=0 



t=0 

w 



dt 



t=0 



$f(0) 



$f , o$fo$- ,/ o^(0) 



$f 4 o $f (0) 



= X«(0)) = X,( ?/ ). 

□ 

In particular, if <5 P = {Id]gn} the Lie group M n equipped with the left invariant sub- 
Riemannian structure associated Xi, . . . ,X m is a Carnot group. In other words, our as- 
sumption implies that any metric tangent cone to (M, d) at p is isometric to a Carnot group. 



30 



Notice that when p is regular it has been shown in [T5] that Q p = {Id^n}. Nevertheless, 
as we see in Example 13. 5\ the metric tangent cone may be a Carnot group even at regular 
points. 

A direct consequence of Proposition 13.41 is that if k p = 2 at every p E M then M is 
equiregular and has a step 2 Carnot group as metric tangent cone at each point. 

Remark 6. Recall that n\{p) is the dimension of span{Xi(0), . . . , X m (0)}. When Q p = 
{IdRn} there exist jx < jl < ■■■ < j m (v) such that Xj 1 . . . , Xj is an orthonormal frame for 
the Carnot group, whereas X^ = for all other indexes k. 

4 The blow-up theorem 

In the next subsections we will always be in the following setup: 

(Al) E is a set of locally finite perimeter in an open set SlcM and p £ JF*E n £l\ 

(A2) o"i, . . . , o m is a local orthonormal frame on O, inducing the vector fields X{ = f o crj, 
<p : Jl — > M n is a system of privileged coordinates centered at p and X\ , . . . , X m are 
defined as in (|4"5l) . 

Note that (A2) is fulfilled by any sub-Riemannian structure (U, /) on M, provided 0, is 
small enough. 

In the previous setup, V, G, d denote the corresponding objects relative to the nilpotent 
approximation (see Section T3.3P and 5 r denote the corresponding dilations. Notice that the 
Lebesgue measure on M. n is well-behaved with respect to the dilations 

6 x (z 1 ,...,z n ) = (X w ^z 1 ,...,X w ^z n ), 

the Jacobian being 

J5 x (z) = \ Q v, Vz e R n 
where (wi(p), . . . , w n (p)) is the vector of weights of V at p (and of V at 0) and 

k p n 

Q p = Y,i&MV l {p)\W-\p)) =$>;(?). 

i=X i=l 

For simplicity, in the sequel we rename Q p by Q and Wi(p) by Wi. 

We are interested in the asymptotic behaviour of 8i/ r ip{E n Q) as r — > 0. Given this 
setup, we can always reduce ourselves to the case when O = W 1 , p = and ip is equal to 
the identity, possibly replacing E by ip(E n 0) and Xi by (see also Proposition I2.2[) . 

This reduction will simplify our notation. In addition, the differential form Cj = (ip" 1 )*^ 
can be written as Cudxx A • • • A dx n with Q smooth and strictly positive in W 1 and we are 
interested in the asymptotic behaviour near the origin. Since the volume form is only used 
to define the divergence, affecting Dx^-e m a multiplicative way (see also the more detailed 
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discussion in the proof of Theorem [23]), we can actually assume that u> = 1 and the measure 
m associated with Co coincides with the Lebesgue measure. 

Before stating our main result, we recall that in a Carnot group G (see Example 13, 1|) 
with Lie algebra g, a Borel set F is called vertical half space if F is invariant along all 
vector fields l£g (i.e., Dx^-f = 0) except a vector field X in V\, for which there is (strict) 
monotonicity, namely Dx^-f is nonnegative and nonzero. Setting v = X(0) £ TqG, we say 
that F is orthogonal to v. 

The following result has been first proved in [24\ Lemma 3.6], see also |14l Proposi- 
tion 5.4] for a different proof (for Carnot groups of arbitrary steps satisfying further alge- 
braic conditions see \33\ Proposition 2.9]). It shows that invariance needs only to be checked 
along directions in the horizontal layer. 

Lemma 4.1. Let G be a Carnot group of step 2, let m be the dimension of its horizontal 
layer V± and let F C G be a Borel set. Assume that V\ contains {m— 1) independent vector 
fields Yi such that DyAp = and a vector field X such that Dx^-f > 0. Then, if Dx^-f is 
not 0, F is a vertical halfspace. 

Recall that k p is the non-holonomic degree of the sub-Riemannian structure at p and Q v 
is defined in (l54|) . 

Theorem 4.2. Under the assumptions in (Al) and (A2) above, the following properties 
hold: 

(a) the family lg 1 .^EnQ) ^ s relatively compact in the L 1 1 oc (]R n ) topology as r — >■ 0/ 

(6) any limit point lp is monotone along the direction X = Yli u E i(p)Xi, i-e~ 

D X 1 F > 0, 
and belongs to the support of D^\p; 

(c) any limit point lp is invariant along all directions X = ^2j_CiXi with q € C°° , 
ET=i c l < 1 and (c(0),v* E (p)) = 0, i.e. 

D X 1 F = 0. 

Moreover, if Q p = {Idign} and k p = 2, then M. n with the left invariant sub-Riemannian 
structure associated with X± , . . . , X m is a Carnot group of step 2 and F is the vertical 
halfspace passing through the origin, normal to ve{p) = l P*f(v Ei (p)). In particular the 
whole family ls 1/rip (Enn) converges to If as r 1 and 

Hm ai4^ = M^ (55) 

no h(B r (p)) JSf"(Si) 
with h(B r (p)) = m(B r (p))/r. 

4 Hereafter, on a Carnot group G we always consider distributional derivatives computed using the 
Lebesgue measure in graded coordinates. 
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Remark 7. Let us mention that when the sub-Riemannian manifold satisfies the condition 

fop = 2 Vp e M, 

then Q p = {Idjjn} at every point and therefore the assumptions in the second part of 
Theorem 14.21 are fulfilled by any finite perimeter set. For instance, this is the case for 
corank 1 distribution (see Example 13. 7J . 

The next remark points out an application of our results to rank-varying distributions. 

Remark 8. Consider the sub-Riemannian manifold of Example 13.51 Outside the plane 
£ = {(x\, X2, x$) | x\ = X2 = 0} the structure is Riemannian, whereas at points x £ £ we 
have k x = 2. Hence, combining the blowup theorem in the Euclidean case with Theorem 14. 21 
above we obtain that any finite perimeter set in this sub-Riemannian manifold admits a 
blowup at each point of its reduced boundary. 

Concerning the proof of Theorem l4.2[ statement (a) is proved in Theorem l4.31 statements 
(b) , (c) are proved in Lemma 14.41 while the second part of Theorem I4.2| which requires in 
addition the assumptions on Q p and k p , is proved in Section [431 



4.1 Compactness 

In this subsection we show that the family of rescaled sets 5i/ r (p(EC){l) is relatively compact 
with respect to the Lj convergence. Here the difficulty in the proof arises from the fact 
that, in some sense, not only the sets but also the metric depends on r, since the rescaled 
sets have finite perimeter with respect to a family of vector fields which does depend on 
r. For this reason, and also because the convergence of vector fields does not occur in 
strong norms, standard compactness results relative to a fixed system of vector fields are 
not applicable. 

Denoting by B\ the closed unit ball relative to d centered at the origin, a simple com- 
pactness argument valid in general metric spaces (see for instance [8]) provides for any 
rj > a partition of B\ in finitely many Borel sets A\ , . . . , A 7 ^,-, and points z^, . . . , z 1 ^^ 
satisfying 

{x | d(x, zV) < |} c AV c {x | d(x, zV)<^} i = 1, . . . , N(r,). (56) 

The proof of the next result is based on the following compactness criterion. Assume that 
Gh C W 1 are Borel sets, and that for any R > and e > there exist r/ = rj(R, e), h(R, e) 
and rrii h G [0, 1] satisfying 

R~ Qy y \1 G -m ih \dx < e for h > h(R, e), (57) 

where A\, . . . ,A V N ^ are as in ([561) . Then (Gh) is relatively compact in the L 1 1 oc (lR n ) con- 
vergence. The proof of the criterion is elementary, since for any R > and e > we can 
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choose ij in such a way that the map 

N(r,) 

l Gh ^ ^2 m iih l A v h>h(R,e) 

i=l 

provides a projection on a compact set (since the rriih are finitely many and belong to [0, 1]), 
e-close in L l (B R ) norm. 

Theorem 4.3 (Compactness). Let E C W 1 be a set of finite perimeter in a neighbourhood 
of 0. Then, if 

\\D a l E \\(B r (0)) , x 

L := limsup - 9 " \ - '' < oo, (58) 

the family of sets 8\/ r E is relatively compact in L* oc (M n ) as r 1 0. 

Proof. By a scaling argument it suffices to show that, for any R > and e > 0, there exist 
rj(R, e) > and f(R, e) > such that 

N{v) , 

(Rr)~ Q / (Is - m r j| dx < e for all r G (0, f), (59) 

7Z[ JSrM*) 



i- 



with m r> j E [0, 1] equal to the mean value of lg on the set S Rr (B 3rjRr / 2 (z^)). 

We choose r\ S (0, 1) satisfying the smallness condition 32®~ 1 cLRri < e, where L is the 
constant in (I58p . c is the multiplicative constant in the Poincare inequality (I23p and TV, 
detailed below, depends only on the (local) doubling constant of d relative to the Lebesgue 
measure. Given 77, because of (j56h and (|5ip . we can find f > such that, for all r E (0,f), 
it holds 

B Rvr /4 (S Rr z?) C S Rr {A r >) C 5 3ft?r/2 (,5 Rr z[0 C B 2Rr (0) . (60) 

Let us check the first inclusion (the proof of the other ones is similar). If d(w,5 Rr z 1 j l ) < 
Rrjr/A, then for r < r(5) sufficiently small from (15 1 f) (see also the equivalent formulation in 
Remark [5]) it holds d(w,5 Rr z^) < R5r/3. Hence d(5±/( Rr }W, z* 7 ) < 5/3, so that (f56|) gives 
$i/(Rr) w ^ ^ and then w G (5fi r (A^). Possibly choosing a smaller f, we can also assume 
that the Poincare inequality (|23p holds at all points 5 Rr z l J with radius 3Rnr/2, for r £ (0,f). 

From ([60]) we deduce that any point belongs to at most TV balls B 3Rrrr / 2 (5 Rr z^), with 
N depending only on the doubling constant of d. Indeed, setting by brevity a = Rnr/4, if 
x belongs the balls BQ a (5 Rr z 1 /) for i G J, then all these balls are contained in B\2 a (x); on 
the other hand, this ball contains the balls B a (5ji T z^), i G J, which are pairwise disjoint by 
the first inclusion in ([60]) . Since B\2 a {x) C -Bisa^Rr^) we get 

J? n (B a (8 Rr z?)) > cB 5 ^ n (iW5«r^)) > c D 5 ^"(i? 12a (5;)), 
where cd is the doubling constant Jl so that J has cardinality at most c 5 D , so that N < c 5 D . 



5 The local doubling property of the Lebesgue measure with respect to the distance d in privileged coor- 
dinates is proved in [37] . 
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Now, using the second inclusion in (|60p and the Poincare inequality (|23p (which, by [3 
is known to hold with A = 1 in length spaces), we can estimate the sum in (|59p with 

^ / \l E -m r>i \dx < y^l \l E -m itr \dx 

= 2^ C ~^n D 9 1 E\\( B 3Rvr/2{SRrZi)) 
i=l 

< ?>cRN W \\D g l E \\{B 2Rr (Q)). 
By our choice of rj, we obtain (|59p . □ 

4.2 Invariant and monotone directions 

Lemma 4.4. Let Q £ L°°(IR n ) 6e a weak* limit point in L°°(M. n ) of ls 1/r <p(Er\n) as r tends 
to and let X = ^iCiXi with c\ E C°°(0) and X^i c ? — 1- Then, if X denotes the 
homogeneous component of order — 1 at of (p*X, the following properties hold: 

(i) if (u E (p), c{p)) = 1, then £ is monotone along X; 
(ii) if {v E (p),c{p)) = then £ is invariant along X. 

Finally, the family lsi /r ip(Enn) ^ relatively compact in the L 1 1 oc (IR n ) topology as r — > 0, and 
therefore £ is a characteristic function If. Moreover, if X is as in (i), belongs to the 
support of D X 1 E . 

Proof. First of all we shall perform the preliminary reduction described at the beginning 
of Section HI so that f2 = W 1 , p = 0, ip is equal to the identity and m = Jzf" '. 

Let us start with a preliminary remark. Since Dx^-e = Yli °iDXi J-s = (c, ||-D fl l^||, 
we can add and subtract v E (0) in the scalar product and use the defining property ([39]) of 
points in the reduced boundary to obtain that 

l i J''^''-^ 1 *l<f(°» = (61) 

under the assumption on X made in (i), while 

\D x l E \(B r (0)) 
W \\D g l E \\(B r (0)) ~° {b2) 

under the assumption on X made in (ii). Thanks to property (|40p valid at points in the 
reduced boundary we have also ||-D g l,E||(-B r (0)) x J?? n (B r (0))/r as r J. 0, and the Ball-Box 
inclusions (|53p give ||-D ff l#||(-B r (0)) x r^ -1 as r J, 0. Hence, using once more the Ball-Box 
inclusions also in the numerators of (161 f) and (]62p , we can write them in the more convenient 
form 

limr 1 ~ Q \\\Dgl E \\-D x l E \(B Rr ) = Vi? > 0, (63) 

rj.0 
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]imr 1 - Q \D x lE\(B Rr ) = Vi? > 0. (64) 

Now we have all the ingredients to prove (i). Fix ip G C^W 1 ) nonnegative and let R 
be such that the ball Br contains the support of ip. By the definition of D^(^, we have to 
prove that 

- f (ijdWXdz- [ ((Xi(j)dz>0. (65) 
Let pi — > be such that lg E weak* converges to ( and define 

Yi := pi{5i/ Pi )*X. 

Recalling that Y\ converge to X and div converge to div X uniformly on compact sets of 
W 1 (see Proposition 13. 4p . it will be sufficient to show that 



i— too 



lim / V divYi dz + / Y^ dz < 



s i/Pi E J5 i/Pi E 
Setting tpi(y) = ip(Si/ Pi y) and changing variables, this is equivalent to 

lim p X r Q ( [ & divX dy+ f Xft dy] < 0. 
\Je Je J 

Now we can integrate by parts, and we are left to show that 

iPidD x l E > 0. 



lim p]- Q f 



B Rpi (0) 



This is an immediate consequence of (|63[) . because ipi are nonnegative, uniformly bounded 
and their support is contained in Bn Pi . The proof of (ii) is analogous, and relies on (|64[) . 

The fact that Q = If for some Borel set F follows by Theorem 14,31 which pro- 
vides compactness in the stronger L 1 1 oc (lR n ) topology (finiteness of L in (f58|) follows by 
\\D g l E \\(B r (0)) x r Q ~ 1 ). In order to prove that belongs to the support of D x lp, under 
assumption (i) on X, we notice that the same argument used above (with integration by 
parts to justify the first equality) gives 

[ X (z)dD x l F (z) = lim f x(z)dD Yz l 5 E (z) (66) 
Jb r ^°°Jbr 1 

= lim p\~ Q \ x(8i/ Pi y)dD x lE(y) 
Jb RPz 

for any x G C^(Br). If we use ([63]) . ||JD fl l£;||(J3 r .(0)) x and assume that x is nonneg- 
ative and x = 1 in a neighbourhood of 0, we get f x(z) D^If > 0, proving that belongs 
to the support of D^lp. □ 
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4.3 Characterization of F when the tangent cone is a Carnot group 

Proof of the second part of Theorem 14.21 Assume Q v = {idiRn}. Proposition 
ensures that ffi n (with the operation *) is a Lie group such that X\, . . . , X m are left-invariant. 
Since the Lie algebra \Ae z {X\, . . . ,X m } is stratified (see Proposition ^, 4p . the group M. n with 
the left-invariant sub-Riemannian structure associated with X±, . . . , X m is a Carnot group. 
Recall that n\{p) = dim'D(O) < m. Define the left invariant vector field Y\ by 

m 

t=i 

Since 1^(0)1 = Gq(ve(0)) = 1, we have G{Y\) = 1. By construction, thanks to Lemma H31 
Y\ is a monotone direction, i.e., D^\p > 0. Let Y2, ■ ■ ■ ,Y ni be left invariant vector fields 

on M. n such that Yi, . . . , Y ni is an orthonormal frame for the Carnot group. Then, again by 
Lemma 14.41 Do lp = for every j = 2, . . . ,n\. Therefore, applying Lemma 14. II we obtain 
that F is the halfspace orthogonal to the geometric normal ^e(O) = /(^(O)). 

Finally, we prove (|55p . After our reduction to the case Vt = R n and ui = dx\ A ... A dx n 
(see Proposition I2.2|) . both total variations ||Z) fl l^;||, ||Z)^l_p|| are computed using Jz? n as 
reference measure. Moreover, thanks to (|5ip . in the left hand side of (|55p we can replace 
B r (p) by B r . Set 

m m 

Yi=j2 u h(°) x " Y 3 = Y. c ^ j = 2 > •••' ni 

i=i 1=1 

where Cjj are such that 1} = X^i=i c ijXi- Thanks to ([oTj) . it holds 



^0 h{B r ) rio h(B r ) 

and, similarly, ||J5^1^ || = Dy^lp. Thus, (also taking (f5Tj) into account) ([55]) is equivalent 
to 

lhn D Yl lE(B r ) = 

By scaling, we can read the property we want to prove as 

limr^V^nlEOBiW) = %1f(£i(0)). (67) 

Now, by (|66p the family of nonnegative measures r 1 ^ (di/^^Dy^E weakly converges to 
D^If as r I 0. Since 9i?i(0) is -D^lir-negligible, applying a well-known convergence 
criterion (see for instance [11] Proposition 1.62(b)]) we obtain (|67p . □ 
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